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I. INTRODUCTION 

For more than two decades there have been studies 
of the nucleation of bubbles during the QCD phase 
transition (QCDPT), the transition from a universe 
of deconfined gluons, quarks and antiquarks (the 
quark/gluon plasma) to a hadronic universe. These 
studies, based on a first-order QCDPT, use mod¬ 
els of QCD to estimate the critical size of bubbles, 
the nucleation rate, and the distance between nu¬ 
cleating centers P, |3,|l,0,[EIEIEIliS ITl| . 
In an attempt to find possible observational ef¬ 
fects from this early universe phase transition, there 
have also been estimates of magnetic fields gen¬ 
erated by nucleating bubble surfaces during the 
QCDPt[ 30) and possible observable CMBR ef¬ 
fects arising from magnetic walls which arise 
from collisions of nucleating bubbles [l5j|. Esti¬ 
mates of magnetic fields arising from bubble col¬ 
lisions during the earlier electroweak phase transi¬ 
tion (EWPT) usin g eff ective Lagrangians also have 
been made fM IitI . Many of these studies make 
use of the Coleman model of tunneling from a 
false to a true vacuum. Nucleation rates were es¬ 
timated using standard thermodynamics/statistical 
mechanics 0 and using hydrodynamics!^ with the 
Langer formalism|^. 

The knowledge that the universe evolved from a 
quark/gluon plasma to our hadronic universe in¬ 
spired the RHIC (relativistic heavy ion collisions) 
program at BNL and other laboratories. The chal¬ 
lenge for the RHIC program, as well as the early 
universe studies, is to identify unambiguous observ¬ 
ables for the transition (or transitions). In a recent 
review [ 2 ^ possible relations between the early uni¬ 
verse and heavy ion heavy phase transitions and the 


status of theoretical attempts to reach thermaliza- 
tion on the time scale that seems to be needed at 
RHIC are discussed. There is also a detailed discus¬ 
sion of flow that is associated with RHIC. In this 
review, as well as in another recent report [^, the 
possible existence of color superconducting phases 
and a color glass condensate that involves low mo¬ 
mentum fields with long time scales compared to 
what is considered to be natural time scales are ex¬ 
amined. Extensive references to theoretical and ex¬ 
perimental research in RHIC physics, and possible 
relationships to the early universe chiral phase tran¬ 
sition, are given in these reviews. 

Lattice gauge calculations indicate that there is 
no first order EWPT transition nor consistency with 
baryogenesis in the Standard EW Model. These cal¬ 
culations also show that with supersymmetric fields 
there can be a first order transition and consistency 
with baryogenesis during the EWPT . Lattice calcu¬ 
lations are not yet able to prove definitely whether 
there is a first order QCDPT. Because of the great 
importance of possible observational effects of a first 
order QCDPT, with bubble nucleation, collisions, 
and possible observable magnetic and other effects, 
in the present note we assume that there is a first 
order chiral phase transition and investigate QCD 
bubble nucleation. 

The most notable aspect of the present work is 
that we start from the basic QCD Lagrangian, rather 
than an effective model, and use previous research 
on nonperturbative QCD condensates to carry out 
the calculations. 

In Sec. H we derive the critical radius for nu¬ 
cleating bubbles, and in Sec. HI we estimate the 
site separation for nucleation of bubbles during the 
QCDPT. 
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II. CRITICAL RADIUS DURING THE 
QCDPT 

In this section we estimate the critical radius, i?c, 
for bubble nucleation, starting from the basic QCD 
Lagrangian. Physically, the critical radius is at¬ 
tained when force due to the pressure difference in¬ 
side and outside the bubble, a volume effect, equals 
that of the surface tension, an area effect. In the 
following section we derive the nucleation probabil¬ 
ity and site separation, which requires quite different 
methods. First, an outline of the paper is given 

A. Outline of Paper 

The starting point of our work is the basic QCD 
Lagrangian density: 

1. QCD Lagrangian Density and basic Eguations of 
Motion 

The QCD Lagrangian density for massless quarks 
is 

£0^^ = -\tt\G ■ G] + qn^{id, + gA,)qf{l) 
where 

G = d^Ai, ( 2 ) 

= A“A 72 , 

with A“ the eight SU(3) Gell-Mann matrices, 
[A“,A^] = 2z/“^°A'^. Minimizing the action, one ob¬ 
tains the general QCD equations of motion (EOM) 

d^d^^Al - d^d^A^ + gr’^%2Ald>^Al 
-A^d^A^'^ - df^A^Al) + 

AlA>^^Al+gqn,^qf = 0 . (3) 

As we shall see, the last term in the EOM, Eq. 0, 
which gives the quark-gluon interaction, plays a cru¬ 
cial role in determining the critical radius of nucleat¬ 
ing bubbles during the QCD chiral phase transition. 


2. Method for Deriving the Critical Radius 

As mentioned above, the force at a bubble wall 
from the pressure difference on the two sides of the 


wall, which drives the nucleation, increases for given 
pressure difference as one extra power of the bub¬ 
ble radius than the force from the surface tension, 
which tends to shrink the bubble. For the calcula¬ 
tion of the surface tension, u, and pressure differ¬ 
ence, Ap, we start with the basic QCD Lagrangian. 
To determine the critical radius we minimize the ac¬ 
tion with res pec t to the radius, as in the models of 
Refs.[lH IliA 1^1. These calculation are done near 
the critical temperature, and do not require finite 
temperature field theory. The approximations used 
are dropping the higher-dimensional gluonic terms 
in calculating the surface tension and using only pi¬ 
ous as hadrons for the effect of the quark condensate 
forming at the critical temperature. These approxi¬ 
mations are discussed in the next subsections in this 
section. 

For this investigation only the action in the vicin¬ 
ity of the bubble wall at a temperature near the 
critical temperature is needed, and we use an exten¬ 
sion of the instanton model to SU(3) symmetry in 
Minkowski space. For the more general study of nu¬ 
cleation during the QCD chiral phase transition this 
model must be generalized to treat the interior and 
exterior of the bubbles. 


3. Method for determining the Nucleation Probability 
and Site Separation 

For the determination of the nucleation site sepa¬ 
ration, the main goal of the present paper, standard 
classical statistical mechanics is used to determine 
the rate at which nuclei of bubbles of critical radius 
form. For this the pressure difference between the 
quark/gluon and hadronic phases at the critical tem¬ 
perature, as well as the surface tension and critical 
bubble size, are needed, and are obtained from the 
QCD Lagrangian density within our instanton-like 
model. 

For the final calculation of the site separation 
distance the temperature dependence of the pres¬ 
sure difference is needed. For this we use the 
standard finite temperature field theory, with a re¬ 
duction from four-dimensional space-time to three- 
dimensional space at fixed temperature by replacing 
real time by imaginary 1/temperature. 

B. SU(3) Ansatz for the color gauge field 

We use the Lorentz gauge and an SU(3) ansatz, in 
which we remove the color dependence of Eq® by 
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extending the color fields, to SU(3) matrices: 
d^Al = 0 

A-l ^ . (4) 

Note that this ansatz cannot be used in the La- 
grangian, Eqs. 12 El), to define the action and obtain 
the EOM Eq.@, but it is rather a method to apply a 
reasonable symmetry on the color field in the EOM. 
It can be compared to the instanton ansatz |^, with 
the color fields given in terms of the SU(2) quantities 
1^ - but our prescription defined in Minkowski 
space leads to Lorentz covariant equations needed 
for studying nucleation and other time-dependent 
processes. We are guided in our choice of parameters 
by the instanton liquid model reviewed in Refill. 
In using this ansatz one operates on Eq. m by the 
color SU(3) matrix A“, uses the matrix extension of 
the color field given by Eq. 0 and the properties of 
the SU(3) generators to obtain the EOM for the 

+ 

-x^ <m> NW, = 0 , (5) 

where we make use of the studyHI of the vector vac¬ 
uum susceptibility defined by the three-point func¬ 
tion, 

< gqiu)A/2q >= x^N <qq>Al, (6) 

in terms of the quark condensate, < qq >. The 
parameter x^ was determined from the study of the 
vector three-point function for a nucleon |27|. while 
the factor N, the number of hadrons in a typical 
hadronic volume V, is needed for our finite T study 
and will be discussed below. 

To determine the critical radius one needs to bal¬ 
ance the forces near the surface of the bubble. The 
essential variable for determining the critical radius 
is s = y/x>^Xfj,, and we use the form 

= x^Wis) , (7) 

with the gauge condition 

W'{s) = --W{s) . ( 8 ) 

s 

The resulting EOM for the function tE(s) is 

W" + ^IT' - ^giw^ + sWW') 

+ ^g‘^s^W^ - x^N <qq>W = 0 . (9) 


C. Surface tension, pressure, and critical 
radins 


Since the purpose of the present paper is to es¬ 
timate the critical radius of nucleating bubbles and 
the site separation, and we are not attempting to 
study the general problem of bubble nucleation and 
collisions during the QCDPT, we shall concentrate 
on the bubble wall and the regions in the QGP and 
HP close to the bubble walls. We neglect some of 
the hadronic structure in the HP and quark/gluon 
structure in the QGP. This allows us to directly use 
the QGD Lagrangian for the nucleation properties in 
the early stages of the QGDPT, which has not been 
done previously, but the method must be extended 
for the complete treatment of the QCDPT. 

The instanton representation of the color gauge 
field is known to give a satisfactory representation 
of mid-range nonperturbative QCD. Although in- 
stantons cannot be defined between hadronic and 
quark/gluon vacua, we shall use an instanton-like 
representation of the bubble wall, as in our earlier 
work [T2[l| . This will lead to a modification of the 
variables used in the expression of the color gauge 
fields, as explained below. 

Let us first estimate the surface tension, one of the 
two parameters needed to determine the critical ra¬ 
dius. The surface tension is obtained by integrating 
the energy density, T°°, through the bubble wall. 


^inst 


/ 


dsT°° 


~ / ds[6W^ - ^g{s^W^ + fs^W^)].(Il) 
7 2 8 


Only the dominant gluonic contribution in Eq. m 
will be used in estimating the surface tension. 

In Refill it was shown that an instanton model 
for the bubble wall is consistent with the surface ten¬ 
sion estimated in lattice QGD calculations. Based on 
this we use an instanton-like form near the surface 
of the bubble wall. Recognizing that Wf^W^ = 0 
at s = So, the position of the bubble wall in the 
instanton model, and is nonzero in the region of ap¬ 
proximately s = so±p, we modify the variables used 
in Eq. 0 ; and for and W{s) we assume the 

instanton-like form 


W{s) 


s^W{sf 

Cw 

(s2+p2)2 ’ 


( 12 ) 
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FIG. 1: Solution for W(s) for t=0 near the surface at 
So = r ~ i?=10.2 fm 


which satisfies the gauge condition Eq. © for s > p. 
Note that with the metric obtained with the con¬ 
tinuation from Euclidean space, in which instantons 
are derived, to Minkowski space, s = — sq; 

the most important region for is the inter¬ 

val —p < s < p. To check the validity of the this 
form we solve the EOM Eq. 0 for W(s), as shown 
in Fig.l. In Fig.l sq = 10.2 fm, which can be inter¬ 
preted as R, the radius of the bubble if we take t=0 
as the time shown in the figure. One can see from 
the figure that the surface peak drops about a factor 
of two as s-R goes from 0 to 0.15, as expected from 
the form of Eq. 113. Therefore, the instanton-like 
form is seen to be valid near the surface of the bub¬ 
ble. For our calculation, we take p = 0.2 fm, which 
is consistent with the result of Fig. 1 and the value 
of the p for the instanton liquid model[^. One can 
find the value of the parameter Cw from such solu¬ 
tions, but as is shown below this is not necessary in 
the present work, where we only derive the critical 
radius and the nucleation site separation. The os¬ 
cillations in W{s) found outside the wall would be 
important for collisions, which we do not treat in the 
present work. See Ref. 

Recognizing that the lowest-dimensional term is 
dominant, the surface tension is 

cr ~ y ds&W{sf = (13) 

To determine the pressure difference, we exploit the 
results of the Higgs model 0, ^|, noting that the 
expectation value of the Higgs field is replaced in our 
work by the quark condensate, which also goes from 
zero to a finite value during the phase transition. In 
the Higgs model, in the thin wall approximation^^ 


m, 

T™(H=r;)-r™(|<).|=0) = eAr?4 (14) 

which is the negative of the last term in the La- 
grangian density used in Refill. Note, however, 
that 

= Ap (Higgs Model). (15) 

Thus, in the Higgs model, as well as for our La- 
grangian, Ap is given by the difference in across 
the bubble wall. With this observation, we are able 
to use all the familiar results of finite-temperature 
field theory for obtaining the consequences of our 
Lagrangian density for nucleation. 

Note that at this stage of the nucleation process 
we are working in the brief time interval between tc 
and tf, the time following tc when the universe has 
reheated (see the following section on nucleation). 
Since the gluonic condensate does not vanish rapidly 
with T ~ Tc, the difference in the gluonic free en¬ 
ergy on the two sides of the bubble wall is negligi¬ 
ble during the first stages of nucleation of interest 
here. For similar reasons we also neglect the free en¬ 
ergy difference of the noninteracting fermionic terms 
within and outside of the bubble, which play an im¬ 
portant part of the QCDPT after the time tf when 
there is an equilibrium nucleation process from the 
quark/gluon to the hadronic universe. For example, 
in a more general treatment of nucleation and colli¬ 
sions for the QCDPT, the perturbative hadronic and 
quark kinetic energy terms in the QCD Lagrangian 
given by qf{i^^d^)qf must be considered |1,|3. 

Since the quark condensate is the parameter for 
the chiral QCD phase transition being treated in 
the present work, as one can see from the expres¬ 
sion for T°°, given in Eq. m, in our theory start¬ 
ing with the QCD Lagrangian the pressure differ¬ 
ence is given approximately by the term in 
that is the negative of the quark-gluon interaction, 
£q(mt) _ _ ^ 5 g 7 ^(A“/ 2 )g > the last term in 
the Lagrangian density O- Having determined the 
surface tension and Ap from our Lagrangian, the 
critical radius can be obtained by minimizing the 
four-dimensional action 

Se = 27rVR3-k y ApR^ (16) 

with respect to the bubble radius, R. 

In applying this formalism to our model we use 
the following concepts: 

i) Our bubble wall is instanton-like, not in¬ 
stantons. This model follows from the success in 
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obtaining an instanton-like interior wall after bub¬ 
ble collisions with QCD bubbles having such an 
instanton-like bubble wall[T5|. 

ii) The wall is very thin, so a thin-wall ap¬ 
proximation is justified. We just need the difference 
between inside and outside the bubble wall. 

iii) Making use of the fact that the quark 
condensate vanishes for T > outside the bubble 
wall and that during the phase transition when T ~ 
Tc inside the bubble pions dominate the hadronic 
phase, we assume that all of the nonvanishing quark 
condensate is in the pions inside the bubble. Us¬ 
ing Eq. ®, the quark-gluon interaction Lagrangian 
density inside the bubble is 

^ -x^n^ <qq> , (17) 

while outside the bubble 

£9(mt) ^ 0, (18) 


would give a 12% decrease in the critical radius and a 
larger nucleation site separation. If one takes a = 2p, 
which is larger than expected for the instanton-like 

wall. Ip = X 0.82. Therefore we conclude that 

the estimate in Eq. m is correct to about 10%. 

Recognizing that the bubble wall thickness is 
given by p = 0.2 /to, we use the thin-wall approx¬ 
imation in minimizing the action 1^ to get the 
classical thin-wall equation for the critical radius. 
From Eqs. il:il2llRH we obtain 


Rc = 


3a 

Ap 


IStt 1 1 

P < 99 > /2 


( 22 ) 


Using the standard value for the pion volume 
densityl^, Utt/U = 0.365T^, with V = (47r/3)i?^ 
= the nucleon volume (i?„=l.l fm), and using the 
value of the vector susceptibility from R ef. 

< qq > /2 ~ 3GeV^, with T — 150MeU, we 
obtain from Eq. 123 


where the color field is not of the instanton-like 
form given by Eqs. 0 ITinil except at the bubble 
wall. A more general theory of gluonic structure in 
the hadronic phase is needed to discuss the interior, 
but for the purposes of the present paper we only 
need the pressure difference on the two sides of our 
bubble wall. 

iv) From this we conclude that the difference 
in inside and outside the bubble wall is 

A[T°°] ~ ^n^<qq><ApA^^> , (19) 

with < ApA^ > the mean value of A^ in the wall. 
Using Eqs. mm, this gives for the pressure dif¬ 
ference inside and outside the bubble surface in the 
thin wall approximation 


Ap 


2 




-n-j, 


< qq> s^W{sf 


< 99 > 


16 p® 


( 20 ) 

( 21 ) 


where we have taken N = n-n, the pion number, and 
have evaluated at s = p. The fact that the pions 
dominate the hadronic density at the phase transi¬ 
tion, and thus that the effective density of quark con¬ 
densate is given by the pion density, is well-known. 

To estimate the error in evaluating (s)^ at s = 
p, we can calculate the mean value by calculating 
la = J_^s'^W{s)‘^/2a. The most reasonable value 

for a is a = p. One finds Ip = x 1.12. This 


Ra = 11/to. (23) 

There are a number of approximations which 
could change our value for and due to the sensi¬ 
tivity of the nucleation site separation to this param¬ 
eter, this could change the evaluation in the next sec¬ 
tion. We show that this value is consistent with the 
evaluation of Ref[3 at the critical time of freezeout 
in the next section. The innovation of the present 
work is that we estimate Rc directly from the QCD 
Larangian, rather than using models as in previous 
estimates. 


III. NUCLEATION PROBABILITY AND 
NUCLEATION SITE SEPARATION 

For applications to astrophysical observations the 
distance between nucleation sites, dn, is a critical 
parameter, since with large separation between nu¬ 
cleation sites the hadronic universe can be formed 
via the collision between a few large nucleating bub¬ 
bles, and this could lead to interesting large-scale 
structure. E.g., it was shown with such a scenario 
that lar ge m agnetic walls would form during the 
QCDPT[T^, which could lead to observable effects 
in Cosmic Microwave Background polarization. 

One obtains the nucleation site separation, d„, 
from the probability of nucleation per volume-time, 
p(t). In this section we estimate the nucleation prob¬ 
ability using the critical radius derived in the previ¬ 
ous section. In classical statistical theory the rate at 







6 


which nuclei form per unit volume-time is given by 

P{Tit)) = , (24) 

where in a Coleman-type model of tunneling from 
the false to the true vacuum the action, S{t), is 
treated in Euclidean space, and for a nucleus of 
critical radius in four dimensions is 5" = 2'K'^aR^ — 
^ApR^, the extremum of which gives the relation¬ 
ship of Eq. 1221). Most of the work that has been 
done during the past two decades has assumed ho¬ 
mogeneous nucleation, for which the form 

P{t) = P{tf)e-°‘^*f-*^ , (25) 

with tf the time at which the universe has reheated 
after the time tc, at which time T is the critical tem¬ 
perature Tc for the phase transition (assumed first 
order). The parameters of the theory are determined 
from the energy momentum tensor in a held theory 
or from the thermodynamic potential in a theory us¬ 
ing classical statistical mechanics. In the following 
subsection we discuss the application of our QCD 
approach assuming homogeneous nucleation, and in 
the next subsection briefly consider inhomogeneous 
nucleation. 


A. Homogeneous nucleation 

From statistical mechanics one knows Q that the 
probability for a huctuation producing a bubble of 
radius Rc for T < Tc satisHes 

P{T) oc , (26) 

with AF the bulk free energy difference between the 
two phases plus the contribution from the surface 
free energy, 

AlTT 

AF = —RliPg - Ph) + 47rai?2 , (27) 

with {Pq,P}i) the pressure in the (quark,hadronic) 
phases, and in our previous notation Ap = Pq — Ph- 
In the estimate of nucleation in the bag model 0, 
0 the quark pressure at T < Tc, with B the bag 
constant, is 

P„ = ^^T^ - B . (28) 

The main difference between our present work and 
that of the earlier work is that we use the QCD La- 
grangian directly and the vanishing of the quark con¬ 
densate at T = Tc, rather than the bag model to 
obtain AF. 


In the approach of 0 Csernai and Kapusta, 
who use the Langer model[^, and take Ph = 
(37r^/90)r'*, find a singularity in Rc{T) as T = Tc, 
and obtain a value of Rc in agreement with our value 
of Rc — 11 fm for {Tc — T)/Tc — 1%. Noting that 
{Tc — Tf)/Tc is of the order of 1%, one sees from 
Eq.||22Jl that our result for the critical radius Rc is 
not inconsistent with this work. 

We now continue with the standard statistical me¬ 
chanics approach^. To obtain the nucleation site 
number density one carries out the integral 


= / p{T{t))f{t)dt 


rB 


Pirn) 


(29) 


where f{t) is the fraction of the universe which has 
not nucleated at time t, and EoCT uses the result of 
Ref|3 that f{t) is a step function quickly disappear¬ 
ing at the time tj shortly after tc- In this picture 
one obtains for the site separation distance 


dn — 0.3 


tV^l 


(30) 


The latent heat density is obtained from the free 
energy difference between the phases byi 

L = Tc-^Ap. (31) 

Since the T dependence of Ap is needed, one can¬ 
not use the result of Eq(|2Il), but must return to 
the thin-wall relationship between Ap and IT^IT^, 
EoCT. Using our observation that in the vicinity 
of the bubble walls the solution for lU^IU^ has the 
instanton-like form of Ealtnil. we find that 


^(^2 +^ 2 ) 2 ’ (32) 

with K a constant. In real time Euclidean space 
s§ = {r — R)'^ + Making the reduction from 4- 
Euclidean space to equilibrium, —IjT'^, and 

s§ = {t — RY — (1/T)^. Using this form and the fact 
that for s\ = , K = 4p^AF — we find that 


L 


4 3ct 
Tf^llc ■ 


(33) 


Using the values p = .2/to, the lattice gauge value 
for a, and our result that Rc = 11.0/to, we find from 


Ea dSni) 


dn = 5.23m. (34) 


As anticipated, this value is about an order of mag¬ 
nitude larger than previous bag model results. 
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B. Inhomogeneous nucleation 

The idea that impurities can have an important 
effect on the probability of nucleation for cosmolog¬ 
ical phase transitions has been known for over two 
decades 1^. Recently this has been considered for 
the QCDPT[^[ii|. As pointed out in , impurities 
in the universe, with number density riin at a time 
ti, tc < ti < tf, can be represented by a P{T{t)) 
which differs from the homogeneous value of Eq. (EH 

by 

P{t) = + . (35) 

Using models for this form the nucleation site sep¬ 
aration can be orders of magnitude larger than the 
value for homogeneous nucleation. 

We would of course obtain similar model- 
dependent results, but since this does not follow di¬ 
rectly from our QCD picture we do not consider it 
further here. 


IV. CONCLUSIONS 

In this work we have started from the energy- 
momentum tensor derived from the basic QCD La- 
grangian with quark and gluonic color fields and 
assumed a first-order QCD phase transition when 
the temperature of the universe is ~ ISOMeU. 
From this the surface tension of gluonic walls in 
an instanton-like SU(3) treatment of the color field 
is derived. The free energy difference between the 
two phases on the opposite sides of the bubble wall 
arises mainly from the vanishing of the quark chiral 
condensate in the quark/gluon phase, and we derive 
this AU using previous work on the nonperturba- 
tive vector quark three-point function. Recognizing 


that the wall is only 0.2 fm thick, the classical action 
is minimized in the thin-wall approximation to de¬ 
rive the critical radius for nucleation, Rc, for which 
the pressure difference driving nucleation overcomes 
the surface tension contracting the bubble. We find 
that Rc — 11 fm, which is satisfactory for pos¬ 
sible astrophysical observations following from the 
QCDPT nucleation and bubble collisions. Assuming 
standard homogeneous nucleation, using our results 
for the critical radius, and deriving the latent heat 
difference density from the instanton model of QCD 
and the quark condensate free energy difference, we 
find that the separation of nucleation sites is of the 
order of several meters, more than an order of mag¬ 
nitude larger than the standard QCD Bag model 
that has been used by previous authors. As research 
by other investigators has shown, the possibility of 
inhomogeneous nucleation could lead to quite large 
distances between nucleation sites and a QCDPT 
with only a few active bubbles. In our next research 
on this topic we shall investigate possible observable 
magnetic structures that would arise from QCD nu¬ 
cleation processes. 
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